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Abstract 

We show that minimal boson stars, i.e. boson stars made out of scalar fields without 
self-interaction, are always classically unstable in 5 space-time dimensions. This is true 
for the non-rotating as well as rotating case with two equal angular momenta and in both 
Einstein and Gauss-Bonnet gravity, respectively, and contrasts with the 4-dimensional 
case, where classically stable minimal boson stars exist. We also discuss the appearance 
of ergoregions for rotating boson stars with two equal angular momenta. While rotating 
black holes typically possess an ergoregion, rotating compact objects without horizons 
such as boson stars have ergoregions only in a limited range of the parameter space. In 
this paper, we show for which values of the parameters these ergoregions appear and 
compare this with the case of standard Einstein gravity. We also point out that the 
interplay between Gauss-Bonnet gravity and rotation puts constraints on the behaviour 
of the space-time close to the rotation axis. 


PACS Numbers: 04.70.-s, 04.50.Gh, 11.25.Tq, 04.20.Jb, 04.40.Nr 


1 Introduction 

Non-topological solitons are distinct from topological solitons pQ in the sense that while the 
latter possess a charge that is of topological origin, the former possess a conserved Noether 
charge that arises from a continuous symmetry inherent in the model. The best known example 
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of a non-topological soliton is the Q-ball, where the Q refers to the conserved Noether charge 
121EUD- The model contains a complex scalar held and possesses a global U(l) symmetry. 
Localized objects in the soliton-sense are possible if a self-interaction potential for the scalar 
held is introdnced. This allows for a snbtle interplay between qnantnm mechanical principles 
and the scalar held self-interaction. In uum non-rotating and rotating Q-balls in (3 -|- 1) 
space-time dimensions have been constrncted nsing a non-renormalizable scalar held potential of 
6th order in the scalar held, while in snpersymmetric extensions of the Standard Model Q-balls 
also exist for more complicated scalar held potentials [SlElDro]. These have been discnssed in 
detail for a scalar held potential of exponential form arising in gange-mediated snpersymmetry 
breaking PSOlE]. 

The self-gravitating connterparts of Q-balls, so-called boson stars have also been discussed 
extensively [121 [131 IHl ESI ESI EZl EEl EH]. In this case, it is sufficient to introduce a mass term 
for the scalar held, while self-interaction of the scalar held is not necessary for boson stars to 
exist. Following the literature, we will refer to these boson stars made out of a massive scalar 
held without self-interaction as to minimal boson stars. These have been studied for the hrst 
time in In most studies of non-rotating and rotating boson stars, however, solutions with 
a hat space-time limit have been discussed, such that the scalar held of the boson star is always 
self-interacting. In 4 space-time dimensions, these solutions have been studied in great detail in 
pm 121] using a 6th order scalar held potential and an exponential potential in m, respectively. 

While boson stars in 4 space-time dimensions are interesting from an astrophysical perspec¬ 
tive ra as well as when considering the discovery of a fundamental scalar held in nature [23] . 
they can also be used to gain more insight into the fundamental properties of the gravitational 
interaction. Most current suggestions for a Quantum Theory of gravity require the existence 
of extra dimensions. Now, one would additionally expect Quantum gravity ehects to become 
important for strong gravitational helds. The gravitational helds of boson stars can be strong, 
but the space-time does not possess horizons. It is hence interesting to study these strongly 
self-gravitating systems and compare their properties with those of black holes. Non-rotating 
boson stars in d-dimensional Einstein gravity with d = 3, 5, 6, 7 have been studied in [2^ using 
an exponential scalar held potential and it was noticed that the critical behaviour of the so¬ 
lutions depends on the number of space-time dimensions. Moreover, non-rotating boson stars 
in 5-dimensional Einstein-Gauss-Bonnet (EGB) gravity have been considered in [25] and it 
was shown that the qualitative features of the solutions change when the Gauss-Bonnet term 
dominates the gravitational interaction. 

The study of boson stars in higher dimensions is also interesting from another point of view: 
while rotating objects in 4 space-time dimensions possess only one angular momentum, they 
can possess more than one angular momentum in more than 4 dimensions due to the existence 
of additional (orthogonal) planes of rotation. In 5 space-time dimensions, e.g., two orthogonal 
planes of rotation exist and hence rotating objects can possess two angular momenta. In the 
case of equal angular momenta the symmetry can be enhanced such that the metric functions 
depend on the radial coordinate only. Rotating boson stars in 5-dimensional space-time with 
two equal angular momenta were discussed for the first time in m using a 6th order scalar 
held potential. It was shown that the sum of the angular momenta is proportional to the 
Noether charge in this case. This study has been extended to include the Gauss-Bonnet (GB) 
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interaction in |26l ED and it was shown that rotating boson stars in EGB gravity do not 
exist when the Gauss-Bonnet interaction dominates the gravitational interaction. Using a 
perturbative expansion it was shown that the solutions cease to exist in this case [2Z]. In 
this paper we point out that this can be traced back to the interplay between the scalar held 
function and the metric functions close to the origin. 

The classical stability of both Q-balls and boson stars is of crucial importance. Gonsidering 
Q to denote the number of scalar particles of which the boson star is made off, the total 
mass M can be compared with the mass of Q scalar bosons with mass m. For M < mQ we 
would expect the boson star to be classically stable in the sense that the kinetic energy of Q 
scalar bosons can be balanced by the gravitational energy of the system. The hrst detailed 
study of this type was done for minimal boson stars in 4 space-time dimensions [28] and it was 
shown that minimal boson stars in 4 space-time dimensions can be classically stable. In all 
studies including self-interaction of the scalar held it was found that stable as well as unstable 
boson star solutions exist (see [2^ for a recent study). This is true for non-rotating as well as 
rotating boson stars in 4 and more dimensions. In general, it was found that the increase of the 
gravitational coupling leads to a decrease in the relative extent of the classically stable branch 
with respect to the classically unstable branch. For sufficiently large gravitational interaction 
the solutions are always classically unstable. 

One of the objectives of this paper is to point out that minimal boson stars in 5 space-time 
dimensions are always classically unstable - both for the non-rotating and the rotating case as 
well as for Einstein gravity and Gauss-Bonnet gravity, respectively. This is remarkable since 
in 4 space-time dimensions stable minimal boson stars exist and thus demonstrates that the 
number of space-time dimensions can inhuence the stability. We point out below that this is 
related to the 1/r^ fall-off of the gravitatioual poteutial as compared to the 1/r fall-off iu 4 
dimeusious. Let us remark that rotatiug minimal boson stars in 5-dimensional Einstein gravity 
have been considered in the context of rotating black holes carrying scalar hair, where they 
appear as limiting solutions [30] . 

Another objective of this paper is the discussion of the ergoregions of the rotating boson 
stars. In general, rotating objects can possess an ergoregion in which the asymptotically time¬ 
like Killing vector becomes space-like. As pointed out iu [SH [32] this cau lead to a superradiaut 
iustability because infalling bosonic waves are amplihed when reflected. For boson stars it has 
been argued that the appearance of ergoregions leads to an instability, the so-called ergoregion 
instability [33]. In contrast to black holes, the space-time is globally regular aud possesses uo 
horizous, such that scattered waves that cau escape to inhnity carry away euergy aud by such 
a process cau destabilize the star. 

In this paper we demonstrate that rotating boson stars in 5 space-time dimensions can also 
possess ergoregions, but only for sufficiently large increase of the scalar held function close to 
the origin. We also show that the GB interaction changes the features of these ergoregious only 
marginally. 

Our paper is orgauized as follows: iu Sectiou 2 we iutroduce the model aud give the Ansatz. 
Iu Sectiou 3, we discuss our numerical results. This iucludes the discussiou of the interplay 
between rotation and the Gauss-Bonnet interaction, of the (in)stability of minimal boson stars 
as well as of the appearance of ergoregions, respectively. Section 4 contains our couclusious aud 
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outlook. 


2 The model 


In this paper we study Einstein-Gauss-Bonnet gravity in 5 space-time dimensions coupled to a 
complex scalar field that can possess self-interaction. The action of this model reads (assuming 
natural units such that h = c = 1) 


S = 




(TX 


R + 


ihTrGs + m^n^n + v;i(n^n)) 




where R is the Ricci scalar, 11 is a complex scalar doublet with mass m and self-interaction 
potential Rsi(n^n), a denotes the Gauss-Bonnet coupling constant and the Lagrangian density 
of Gauss-Bonnet gravity reads 


£ 


GB 


dMNKL d 
— n rtMNKL 


AR^^Rmn 


+ R^ 


( 2 ) 


with M, N,K,L E {0,1, 2, 3,4}. Newton’s constant in 5 dimensions is related to the Planck 
mass Mpi^s and Planck length /pi^s in 5 dimensions, respectively, by G^ = = /pjg, while 

the Gauss-Bonnet coupling a has the dimension of a [length]^. If we would treat the action 
above as a low energy effective action of String Theory, a would fulfill a ~ /pj 5 . 

The equations of motion then read 

d 

Gmn + -^Hmn = ^ttG^Tmn , M, iV = 0 , 1 , 2 , 3,4 , (3) 

where Hmn is given by 


Hmn — 2 [RmabcRn^^ ~ ‘^RmanbR^^ — ‘^RmaRn + RRmn) 

- ]^gMN {R^ - ^RabR^^ + RabcdR^^^’^) , R,5,G = 0,1,2,3,4 , 

and Tmn denotes the energy-momentum tensor of the scalar field 

dC 


Tmn — gMNT — 2 

T 


dgMN 


= + dNU^dMU 


— gMN 


g^^ {dKH^dLU + dLU^dKH) + + Ki(nin) 


(4) 


(5) 


The scalar field equation is given by the Klein-Gordon equation 

= ° ■ («> 

In this paper, we are interested in two cases: (a) boson stars composed of a massive scalar field 
without self-interaction, i.e. Ki = 0 (following the literature we will refer to these boson stars 
in the following as minimal boson stars) and (b) boson stars with a self-interaction of the form 

vi,(|nn = mVE(-i)"‘(y) . P) 
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where r/ is a (dimensionful) coupling constant. This self-interaction potential appears in gauge- 
mediated supersymmetry breaking with breaking scale t] [3101. We will concentrate on case 
(a), but also present some results for case (b). 


2.1 Ansatz 


In principle, localized objects in 5-dimensional space-time can possess two independent angular 
momenta associated to the two orthogonal planes of rotation. If one restricts to the case of 
equal angular momenta the symmetry of the object is enhanced and the Ansatz for the metric 
reads [ 21 ] 

ds^ = —bir)dt^ + -^r^dr'^ + g(r)d9'^ + hir) Tsin^ 9 idgDx — Wir)dt)^ -\- cos^ 9 (d(p 2 — W 
f{r) 

+ {g{r) — h{r)) 9 cos^ 9{d(pi — d(p 2 )‘^ , (8) 


where 6 ^ G [0 : 7 r[, while <^ 9 ^ G [0 : 27r[, k = 1,2. The corresponding space-time possesses two 
rotation planes at 6 * = 0 and 6 * = 7 r /2 and the isometry group is M x 17(2). The metric ([ 8 |) still 
leaves the diffeomorphisms related to the definition of the radial variable r unfixed. This can 
be fixed by choosing g{r) = and we will employ this choice in the following. 

In order to construct rotating boson stars in 5 dimensions, the following Ansatz for the 
complex scalar doublet was first introduced in [ 21 ] and reads 

U{t, r, 9, (pi,(p 2 ) = 0(r)e“^*n(6', <pi, (pa) , (9) 


where 11 is a doublet of unit length that depends on the angular coordinates only and is chosen 
such that 

n = (sin0e*‘^bcos0e*‘^=)‘ . (10) 

Note that the case of non-rotating Gauss-Bonnet boson stars, which has been studied in 
[25] . corresponds to the choice II = (1,0)7 In this case, the equations lead to W{r) = 0 and 
h{r) = r^. 

While the metric (jS]) has three commuting Killing vector fields dt, [21], the scalar 

doublet of the rotating solution with fllOp is only invariant under one possible combination of 
these vectors, namely under dt — co (9<^^ -|- d^^). 


2.2 Physical quantities 


The action possesses a global U{1) symmetry which leads to the existence of a globally conserved 
Noether charge Q which reads m 


Q = 



( 11 ) 


where corrsponds to the f-component of the locally conserved Noether current = 
—i — S^Ifln). inserting the specific Ansatz into flTT]) we hnd in our choice of met¬ 

ric: 


Q = 



{u + W) (j)^dr . 


( 12 ) 
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Using the Komar expressions to evaluate the angular momenta Ji = J 2 = J it was realized 
that Q = \ Ji\ + IJ 2 I = 2\J\ [2i] . 

The mass of the solutions can be given by the relevant Komar expression (see e.g. [21]) and 
is proprtional to the prefactor in the 1/r^ fall-off of the metric function b{r) at infinity (see e.g. 

[SI). 

The Ricci scalar, which describes the local scalar curvature of the space-time reads 


R 


8 2h / h” h'2 

2 r4: \ 7- r> 1 .9 


h 2h‘^ rh 


2h' 4 b'h' 

-^ + 2 - 


2bh 


' 2h r 26 


26' 6" hlU'2 6'2 \ 

l^~~b^ 26 ^2^y 

( 13 ) 


where the prime denotes the derivative with respect to r. 

Recall that in pure Einstein gravity, i.e. for a = 0 we have 3R = — IbvrGsT, where T 
denotes the trace of the energy-momentum tensor. 


3 Results 


In the following, we will discuss our numerical results that we have obtained by solving the 
coupled system of ordinary differential equations by employing a Newton-Raphson adaptive 
grid iteration scheme |35j . 

We apply the following rescalings 

a 

—y —, , M = 0 , 1 , 2 , 3,4 , oi —y — , oj —y ttiuj , (p —y Xcp , ( 14 ) 

m m 

where A is a scale. The system of equations then depends only on a, u and n := SttGsA^ and 
additionally on the rescaled self-coupling in the non-minimal case. We then choose A 

as follows 


1 . for minimal boson stars A 


1 

v^SttG's 


such that K 


1 , 


2. for self-interacting boson stars X = rj such that k = 


3.1 Interplay between rotation and Gauss-Bonnet interaction 

In order to understand how the interplay between rotation and the Gauss-Bonnet term effects 
the properties of boson stars, let us hrst remind the reader of the qualitative differences of 
non-rotating boson stars in Einstein and Gauss-Bonnet gravity, respectively. Non-rotating 
boson stars in 5-dimensional Einstein gravity were hrst studied in HD- These share many 
features with the 4-dimensional counterparts: they exist down to a minimal frequency Wmin > 0 
from where several new branches of (unstable) solutions appear that typically form a spiralling 
behaviour. The critical solution at ojcr > corresponds to a solution that has the central 
value of the scalar held function, 0(0), tend to inhnity. In contrast to that Gauss-Bonnet 
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boson stars do not have a minimal value of u. The spiral present in the former case unfolds 
and solutions exist down to a; = 0 if a is large enough and hence the GB term dominates the 
gravitational interaction. However, the value of 0(0) is now restricted and stays below a critical 
value 0(O)cr < C )0 |25]. This seems to suggest that while for Einstein gravity we can localize 
the scalar field and hence the energy density arbitrarily close to r = 0 , this is impossible for 
the Gauss-Bonnet case. We believe that this is related to the fact that Gauss-Bonnet gravity 
has a fundamental minimal length scale inherent in it, namely the Planck length /pi 5 which is 
related to the parameter a. 

Rotating boson stars in Einstein gravity show an analogue behaviour as compared to the 
non-rotating solutions [211 Eg. In the case of rotating solutions, 0(0) = 0 and hence the 
derivative of 0(r) at zero, 0^0), can be used as a parameter. Again, several branches exist and 
the solutions exist in a limited range of u with cU min > 0 , while solutions can be constructed 
for arbitrarily large values of 0^0). Now, it has been observed in [Sg that this is different for 
the Gauss-Bonnet case. In this latter case the value of 0'(O) is limited to a hnite critical value 
at which the solutions cease to exist as was pointed out in Here we demonstrate that the 
criticality of the solutions appears already on the level of the behaviour of the functions at the 
origin. For that we expand the functions close to r = 0 taking the boundary conditions into 
account: 


/(r -C 1) = 1 + F 2 r‘^ + 0{r‘^) , b{r <C 1) = Bq + B 2 r‘^ + 0{r^) , h{r -C 1) = r‘^{l + H 2 r‘^+ 0{r‘^)) , 

(15) 

W(r « 1) = W„ + + 0(r‘) , ^(r « 1) = ^ir + 0(r®) , (16) 

where F 2 , Rq,- 82 , Wo, H0, 0i are constants to be determined numerically. Note, however, that 
the equations of motion lead to several relations between these constants, namely 


^ ft0i + 3(^2 H 2 ) ^ ^01 (Wq u ) 

' ° 3a(F2 + ifa) - 3 ’ ' 60(^2 + 31/2) - 6 

as well as 

2 F 2 H 2 5 H 2 ) — 6(^2 T H 2 } — 2k,<P\ = 0 . 

This latter equation can be solved for either H 2 and F 2 . The solutions for H 2 is 


(17) 

(18) 




5 5a 15a 


with discriminant A given as follows 


(19) 


A = — 12a^F2^ -I- 24aF2 + lOaK0^ -1- 3 


( 20 ) 


7 






Figure 1: We show the quantity A/3 (solid, see fl2U]) l for different values of 01 in dependence 
on a for rotating minimal boson stars. We also give the corresponding mass M (short-dashed) 
and the charge Q (long-dashed). 


In Fig. [T]we show the value of A/3 in dependence on a for different values of 0i. We also 
give the corresponding mass M and the charge Q of the solution. As indicated by fl20|) we find 
that the larger 0i the smaller is the value of a at which A = 0. This is related to the fact 
that at r = 0 the energy-mass density is dominated by the kinetic term encoded in 0i and 
hence, when a becomes large, Gauss-Bonnet gravity dominates the gravitational interaction. 
And since Gauss-Bonnet gravity has a natural minimal length scale encoded in it, we would 
expect solutions ceasing to exist for sufficiently large 0i. 

From the data in Fig. [T] we find the critical values of a and corresponding values of M, 
Q and uj as given in Table [U The smaller we choose 0i the larger is the value of a at which 
A —)■ 0. Our data suggests (which makes also sense for dimensional reasons) that 


ar 


2 

30 ? 


( 21 ) 
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01 

Clcr 

M 

Q 

OJ 

0.5 

2.36 

140.0 

116.3 

0.93 

1.0 

0.67 

178.1 

155.7 

0.99 

2.0 

0.18 

201.8 

180.7 

0.97 


Table 1: Value of critical a for given value of (j)i and corresponding values of M, Q and u. 



r 


Figure 2: We show the metric functions h{r)/r‘^ (solid) and /(r) (short-dashed) as well as 
the tt-component of the metric, —gu (dotted-dashed) and the scalar held function 0(r) (long- 
dashed), for 01 = 0.5 and a = 0 (red) and a = 2.35 (black), respectively. Note that a = 2.35 is 
close to the critical value of a at which A —)• 0 (see also Fig. [T]and Table [T]for detailed values). 


In Fig. [2] we show the behaviour of the metric functions for 0i = 0.5 and a = 0 as well as 
a = 2.35 ~ ttcr- As is clearly seen from this hgure, the metric functions /(r) and h{r) change 
their behaviour close to r = 0. While for a = 0 the metric function /(r) decreases from unity 
close to the origin, it increases in the case of strong Gauss-Bonnet coupling. The metric function 
h(r)/r^ shows the opposite behaviour: for a = 0 it increases from unity, while in the Gauss- 
Bonnet case it decreases. This can also be seen in Fig. [3l where we show the metric functions 
as well as the scalar held function 0(r) for 0i = 2 and a = 0 as well as a = 0.175 ~ While 
the qualitative behaviour of IF(r), gu and 0(r) close to the origin does not change, the metric 
functions /(r) and h{r)/r‘^ show qualitatively different behaviour. This is a clear indication 
of the fact that for large enough a the Gauss-Bonnet interaction dominates the gravitational 
interaction and hence the scalar curvature is no longer given by the energy-momentum content 
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r 
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Figure 3: We show the metric functions h{r)/r‘^ (short-dashed, upper plot) and /(r) (solid, 
upper plot) as well as the tt-component of the metric, —gu (solid, lower plot), the scalar held 
function 0(r) (short-dashed, lower plot) and the metric function W{r) (solid, lower plot) for 
01 = 2 and a = 0 (red) and a = 0.175 (black), respectively. Note that a = 0.175 is close to 
the critical value of a at which A —)• 0 (see also Fig. [T]and Tabled] for detailed values). 
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of the space-time as in Einstein gravity where R T. This is clearly seen in Fig. |H where 
we compare the local scalar curvature R with —2T/3 for a = 0 and a = 0.175 close to the 
maximal possible value of a, respectively, for 0i = 2. For the former case we know from the 
Einstein equation that R = —2T/3 (letting SttGs = k = 1.0) and we confirmed this equality 
numerically (and hence checked our numerical procedure to be valid). For a = 0.175 the scalar 
curvature R and —2T/3 differ close to the origin. While the local maximum of the trace of the 
energy-momentum tensor, T, is located roughly at the maximum of the scalar field function 
0(r) (see Fig. 0]) and the increase of the Gauss-Bonnet interactions increases the r at which 
the maximal energy-momentum content is located, the scalar curvature i? at r = 0 decreases 
with increasing a and shows no longer a pronounced minimum at some hnite value of r. 



r 


Figure 4: We show the tt-component of the metric, —Qu (dashed), as well as the Ricci scalar 
R (solid) for a = 0 (red) and a = 0.175 (black), which is close to the maximal possible value 
of a for which A —)■ 0 for the chosen value of = 2 (see also Fig. [1]). For a = 0 it follows 
from the Einstein equation that R = — 2T/3 (with our rescalings such that k, = 1), while we 
also show —2T/3 (black dotted-dashed) for a = 0.175 to demonstrate that the local energy- 
momentum no longer determines the local scalar curvature in the case where the Gauss-Bonnet 
term dominates the gravitational interaction. 

3.2 Instability of minimal boson stars in 5 dimensions 

The biggest fraction of the results presented in this subsection correspond to that of a massive 
scalar held without self-interaction. This case - at least to our knowledge - has not be studied 
in detail in the literature so far. As pointed out above, we can rescale coordinates and helds 
such that we can set k = m = 1 in this case, such that the only free parameter is the angular 
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frequency oj in Einstein gravity. For Gauss-Bonnet gravity a is a second free parameter. 



CiJ 



(a) M, Q (b) (/.(O) , 6(0), R(0) 


Figure 5: We show the mass M (black) and the charge Q (red) in dependence on the frequency 
oj G [0.9 : 1] for non-rotating minimal boson stars and for a = 0 and a = 1, respectively (left). 
We also show the dependence of the value of the scalar held function 0(r) at the origin, ()i)(0) 
(solid, red), of the value of the metric function 6(r) at the origin, 6(0) (solid, blue), of the value 
of the scalar curvature R at the origin, i?(0) (dotted-dashed, black), as well as of the mass M 
(solid, black) on a; G ]0 : 0.9] for a = 1 (right). 

3.2.1 Non-rotating solutions 

In contrast to the rotating case, non-rotating solutions have a non-vanishing value of the scalar 
held at the origin, 0(0), while the derivative vanishes there. Hence, non-rotating boson stars 
are typically characterized by the value of the scalar held at the origin, which is determined by 
the value of u and vice versa. 

We hrst discuss the case 0 = 0. We hnd that minimal boson stars exist for u G [0.9477 : 1]. 
This frequency does not uniquely characterize the boson stars. Indeed several solutions exist 
with the same value of cu, however, with diherent values of M and Q. This is shown in Fig. 
0 (left), where we give the mass M and the charge Q as function of u. Several branches of 
solutions are present (as in the case of self-interacting boson stars) and the typical spiraling 
behaviour is seen. 
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(a) M, Q, <j>{0), 6(0) 


(b) 6(0), fm, RiO) 


Figure 6: We show the mass M (solid, black), the charge Q (solid,red), the value of the metric 
function b{r) at the origin, 6(0) (solid, blue), as well as the value of the scalar held function 
(j){r) at the origin, 0(0) (solid, green), in dependence on a for non-rotating minimal boson 
stars with u = 0.9 (left). We also give the value of the metric function b{r) at the origin, 6(0) 
(solid, subhgure), the minimal value of the metric function /(r), fm (red, subhgure), as well as 
the value of the scalar curvature R at the origin, i?(0) (black), in dependence on a for these 
solutions (right). 

The limit 0(0) —)■ 0 corresponds to a; —)■ 1 and the scalar held function approaches the 
null function uniformly. Interestingly, the mass and charge are hnite in this limit = 

Quj^i ~ 226.0. This was also observed for boson stars with self-interaction in |23] and the 
presence of this “mass gap” is a special feature in 5 space-time dimensions. Boson stars exist 
for arbitrarily large values of 0(0), however the mass, the charge and u remain hnite. The 
spiral ends at 0(0) oo with ui ~ 0.96, M ~ 78.95 and Q ~ 73.7. 

The scenario is diherent for the Gauss-Bonnet case {a > 0). In the region a; ~ 1, (i.e. 
for 0(0) <C 1) the solution is rather insensitive to a, which is natural since the local energy- 
momentum content is small. However, quantitative and qualitative diherences occur when 0(0) 
is increased. As shown in Fig. Elfor a = 1 the spiral has disappeared and solutions now exist 
for all a; G ]0 : 1]. We also show some parameters characterizing the solution in Fig. |5] (right). 
As can be seen in this hgure, we hnd that the value of the metric function 6(r) at the origin, 
6(0), tends to zero. This suggests the appearance of a singularity at r = 0 in the limit a; —)■ 0. 
In addition, we observe that the Ricci scalar R becomes very large in this limit, while the mass 
stays hnite. 

In order to understand this critical behaviour better, we have also studied the case of a 
hxed bj and varying a, which, in fact, turned out to be more feasible numerically. Our results 
for uj = 0.9 are shown in Fig. |6l where we show the dependence of the mass, charge and several 
parameters characterizing the solution on a. As is clear from this hgure, the mass and charge 
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decrease when decreasing a for fixed u. At the same time (pi increases. In order to understand 
the critical behaviour we have also plotted the value of the metric function b{r) at the origin, 
6 (0), as well as the minimal value of the metric function /(r), fm- /m —t 0 would indicate the 
formation of an extremal black hole which would carry “scalar hair”. However, as was shown 
in [36] the near-horizon AdS 2 x geometry of an extremal Gauss-Bonnet black hole does not 
support scalar hair. Hence the limiting solution is a singular solution with 6(0) 0. 

From the hgures presented above it is obvious that all the solutions without self-interaction 
have M > Q which indicates that these solutions are always classically unstable. This contrasts 
with the 4-dimensional case, where it was shown that minimal boson stars are classically stable 
[28] . Following the arguments in this latter paper it is easy to show why this should be the 
case. Consider the boson star made out of Q scalar quanta of mass m. These scalar quanta 
have each kinetic energy p ~ R~^, where p is the linear momentum of the scalar 

quantum, A its average wavelength and R the radius of the boson star. The kinetic energy 
of the boson star made out of Q quanta is hence Ek Q/R M/{mR), where we assume 
that M ~ mQ. In 4 dimensions the gravitational energy is i?g ,4 ~ —GiM'^/R, where is the 
4-dimensional Newton’s constant. It is thus possible to hnd an equilibrium between the kinetic 
energy and the gravitational energy. As pointed out in [28], the kinetic energy makes the star 
expand until the gravitational energy dominates the system and - because of its attraction - 
allows to have a stable star. Now, this is different in 5 dimensions, where the gravitational 
energy is F'g^s ~ —G^M’^/B?. At small R the gravitational energy will always dominate and 
hence stable conhgurations are not possible. This changes when one includes a repulsive self¬ 
interaction which can balance the gravitational attraction and allow for stable boson stars. We 
show the ratio M/Q in Fig. [7] for non-rotating boson stars without self-interaction and self¬ 
interaction potential ([7]), respectively. It is clear that without self-interaction the ratio M/Q 
is always larger than unity and hence the boson stars are always classically unstable. We also 
demonstrate that the value of the Gauss-Bonnet coupling does not change anything as far as 
this conclusion is concerned. Now, this changes when the self-interaction is present. For k, small 
stable solutions exist, while for sufficiently large k the solutions become again unstable. This 
latter observation is related to the fact that for k too large the gravitational energy dominates 
the kinetic energy (as argued above) and hence stable conhgurations are not possible. The 
hgure demonstrates also that the inclusion of the Gauss-Bonnet interaction does not change 
much for the self-interacting case either. 
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(j, (0) 


Figure 7: We show the ratio of the mass M and the charge Q, M/Q in dependence on the value 
of the scalar field function 0(r) at the origin, 0(0), for non-rotating boson stars with a = 0 and 
a = 1, respectively. We compare the case of boson stars without self-interaction (solid) and 
that with self-interaction potential 14i (dashed). Note that in the case without self-interaction 
we can scale k = 1 without loss of generality, while we have adapted the re-scaling (j) ^ rjcj) for 
the case with self-interaction and hence study different values of k. 
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Figure 8: We show the mass M (black) and the charge Q (red) in dependence on the frequency 
u for rotating minimal boson stars and for a = 0 (solid), a = 0.5 (long-dashed) and a = 1.0 
(short-dashed), respectively. Note that the angular momentum is related to the Noether charge 
like J = Q/2. 

3.2.2 Rotating solutions 

In Fig. [8] we show the mass M and charge Q for rotating minimal boson stars and different 
values of a. The qualitative pattern is very similar to that of non-rotating solutions. For a = 0 
we find the typical spiraling behaviour, while for large values of a the solutions cease to exist 
at some critical value of 0i, respectively u. 

As far as the classical stability of rotating minimal boson stars the conclusions are very 
similar to the non-rotating case. As can be seen from Fig. [8] we find that independent of the 
choice of a the solutions have M > Q and are hence classically unstable. While we would 
expect a repulsive centrifugal force to help balance the gravitational attraction in this case, 
this is not sufficient to render the solutions stable. 


16 











01 

a 

n 

To 

^rgo 

0.5 

0.00 

n.a. 

n.a. 

n.a. 

0.5 

2.36 

n.a. 

n.a. 

n.a. 

1.0 

0.00 

0.41 

0.92 

6.43 

1.0 

0.67 

0.61 

0.98 

6.48 

2.0 

0.00 

0.09 

0.82 

6.10 

2.0 

0.18 

0.09 

0.86 

6.30 


Table 2: Value of the inner and outer radius of the ergoregion, r* and Vo, respectively as well 
as the proper volume of the ergoregion 14rgo for different values of 0i and a. Note that the 
non-vanishing values of a are close to at which A —)• 0 for the respective value of 0i. The 
abbreviation “n.a.” implies that the space-time does not possess an ergoregion. 


3.3 Ergoregions 



r 


Figure 9: We show the tt-component of the metric for different values of and u, respectively, 
and for minimal boson stars with a = 0. 

Rotating black holes typically possess ergoregions between the static limit surface and the event 
horizon. This is the region in which the asymptotic time-like Killing vector held becomes space- 
like and processes such as the energy extraction from black holes (Penrose process) become 
possible. Rotating objects without event horizon can also possess such ergoregions under specihc 
conditions. 

We hnd that these ergoregions also exist for our solutions. This can be seen in Fig. |2]for 
(pi = 0.5 and Fig. IHfor <pi = 2 , respectively. As is obvious from this latter hgure, the value 
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of a does not influence the ergoregion much. It shifts the radius of the ergosphere to slightly 
larger values of r, but otherwise does not influence the qualitative shape. However, it is the 
value of 01 that triggers the existence of an ergoregion. For 0i = 0.5 no ergoregion exists, while 
for 01 = 1.0 it is present already and persists to be present for 0i = 2. We have investigated 
this in more detail. In Fig. |9]we show the tt component of the metric for different values of 
01 and uj, respectively. We hnd that when increasing 0i an ergoregion appears for sufficiently 
large 0i. For a = 0 this happens on the second branch of solutions. 

In Table [2] we give the value of the inner and outer radius of the ergoregion denoted by 
Vi and To, respectively, for two different values of 0i and for the Einstein limit {a = 0) and 
the maximal possible value of a, ctcr, at which A —)■ 0. For both values of 0i we observe that 
the value of Tq is slightly increased. For 0i = 1.0 the inner radius is pushed outwards close to 
the critical value of a, which is also the maximal possible value. Here, the GB interaction is 
strongest and its influence gets stronger when approaching the origin. For 0i = 2 the maximal 
possible a is relatively small, hence we would expect the GB interaction to play little influence 
here. For a hxed value of a we observe that an increase of 0i leads to two things: (a) a shift 
of the inner and outer radius to smaller values of r and (b) the increase of the extend of the 

To _ 

ergoregion in r. When computing the proper volume of the ergoregion Wrgo = / \/—^d^x, 

ri 

where 5^4 is the determinant of the spatial part of the metric, we notice that this increases 
slightly with a and decreases with increasing 0i. For a = 0 we hnd that the ergoregion hrst 
appears at 0i ~ 0.842 when increasing 0i from zero. The volume Wrgo then shows a sharp 
increase from zero and stays around the value 6 for 0 i > 1 . 0 . 


4 Summary and Outlook 

In this paper we have focused our studies on minimal boson stars in 5 space-time dimensions. 
We have pointed out that these minimal boson stars are always classically unstable - in contrast 
to their 4 dimensional counterparts. This is due to the fact that the gravitational interaction 
always dominates at small r and hence a stable conhguration is not possible. Moreover, the 
boson stars possess ergoregions for sufficiently large increase of the scalar held function at the 
origin. We observe that the Gauss-Bonnet interaction alters the location and appearance of this 
ergoregion only marginally. This is also related to the fact that rotating Gauss-Bonnet boson 
stars exist only up to a maximal value of the Gauss-Bonnet coupling and hence solutions with 
arbitrarily large Gauss-Bonnet interaction do not exist. These ergoregions make the solutions 
suher (additionally) from an ergoregion instability. 

It would also be interesting to check whether minimal boson stars in dimensions higher 
than 5 are classically unstable. We believe they are because the power of the fall-oh of the 
gravitational potential increases with the number of dimensions, but leave this for future work. 

In addition, one can construct boson stars with non-equal angular momenta and check 
whether the solutions possess ergoregions. These ergoregions - if they exist - will be different 
in shape to the ones discussed here. The ergoregions of rotating boson stars with equal an¬ 
gular momenta are spherical hypershells, while those of non-equal angular momenta solutions 
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will have an angular dependence. It would be very interesting to see whether the ergoregion 
instability is present and what effects it has. 
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